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Abstract 

Neutrinoless double beta decay (/9/d)ojy occurs through the magnetic coupling of di¬ 
mension five, \w/^v*i among the excited electron neutrino z/*, electron and W boson if 
V* is a massive Majorana neutrino. 
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where Xw is the relative strength, is the composite scale and is the mass of excited 
neutrino. If the coupling is not small, i.e., > 1 and nii,^ = m*, we hnd m* > 3Amw 

which is the most stringent limit. 
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1 Introduction 


If neutrinos are composite particles, there exist the excited neutrinos which couple to the 
ground state leptons by the dimension hve magnetic coupling[1],[2], This interaction is 
expressed as [3] 

Lint = + r]*^L)v*d^W~ + h.c., ( 1 ) 

rrin* 

where v* is a heavy excited electron neutrino, L = (1 — 75 ) 72 , R = {1 + 75 ) 72 , is the 
mass dimension which is of order the mass of z/*, i.e., m*. This interaction is derived by 
the SU{2) X [/(l)gauge invariance and parameters 77 ^ and rjn are normalized by 

17L r +1 VR r= 1 • (2) 


The extensive search of u* have been made by many groups[4] and found that > 
91GeV by assuming that > 1 which is the coupling for z/* —eZ decay similarly 
dehned to Xw- 

The purpose of this paper is to explore the mass range niu* > mz hy using the neu¬ 
trinoless double beta decay {PP)ou by assuming that u* is a massive Majorana neutrino. 
Then the {PP)ou decay occurs through u* exchange. Since z/* enters as a virtual state, 
we can investigate heavy u*. 

In our previous paper[5], we started from the interaction in Eq.(l) and derived the 
effective four point interaction between leptons and hadrons as 


L,ff = -G,,fea'‘-{r,lR + r,lL)v‘d,Jl + h.c., 


( 3 ) 


where J.t is the hadronic current and 


Geff - 2Gj 


ITlj,* 


( 4 ) 


From this effective interaction with the chirality selection rule = 0, the half-life of 
the neutrinoless double beta decay was derived[5] and is given by 
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which disagreed with the old result by Panella and Srivastava[6]. By comparing this 
formula to the Heidelberg-Moscow data[7] for ^®Ge decay, we found[6] 


d!!!!') I + 7^2. |< 1.4.10- 
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( 6 ) 


The formula in Eq.(5) and also the constraint in Eq.(6) are valid as far as we consider 
the effective interaction in Eq.(3). 

Recently, Panella, Carimalo, Srivastava and Widom pointed out [8] that propagators 
of W-boson exchange should not be contracted as far as one sticks to the interaction in 
Eq.(l). This is true since the mass of excited neutrinos m* is greater than mw- They 
calculated[8] the half-life formula of the neutrinoless double beta decay and claimed that 
my calculation is wrong. Here, I say that my calculation is correct as far as we consider 
the effective interaction in Eq.(3). Anyway, I feel to be obliged to repeat the computation 
starting from the interaction in Eq.(l). 

In Sec. 2, we show the computation of the decay and give the half-life formula. In 
Sec. 3, the numerical analysis will be given. Summary is presented in Sec.4. 


2 Decay formula of neutrinoless double beta decay 

In the fourth order perturbations of the interaction in Eq.(l), the {PP)ou decay takes 
place and the S-matrix for this decay is given by 




^-iq{xi-X2) ^-iqi{xi-X3) ^-iq2{x2-X4) 
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(7) 

where is the charge conjugation of e, i.e., = Ce^. 

In the following computation, we take the S-wave function for electrons which is given 
by 
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where e is the energy of electron and Fq{Z, e) is the relativistic Coulomb factor defined in 
Eq.(3.1.25) in Ref.9. The S-wave function is independent of the space coordinate. After 
integrating xi, X 2 and then qi, q 2 , we hnd 


Sfi = -i- 
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where ki = (e*, 0 ). 

Next we perform the integration. There are six poles at ±(i?* — A), ei ± {Ew — E) 
and —62 ± {Ey/ — ie). The poles in the lower-half plain contribute when Xg — > 0 and 

those in the upper-half plain do when Xg — < 0. By performing the integration, the 

time ordering for hadronic currents is dealt automatically. We hnd 

+ VrL) + (£ 2 ) 

x[6'(a;3 -0:4) < Nf\Jl{x 3 )\n >< n\Jl,{x 4 )\Ni > 
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where 
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„ ^ {q* + 

“ E,[{E, + e,)^-E^][{E,-e2y-E^] 

{qw + ki + A; 2 )^gw"^'e*(^'^+^ 2 )(xg-a; 2 ) 

^Ew[{Ew + €2)2 - i?2][(EH^ + ei + €2)2 - 
qwMw -ki- 

^Ew[{Ew - ei )2 - E^Ew - e, - €2)2 - ^ ^ 

Here E^ = ^ml + g* = (i?*, ^ and g* = (i?*, — ^ and similarly for Ew, qw and qw- 
Then, we perform the and x^ integration. The X 3 and X 4 integration can be made 
immediately since we use the non-relativistic approximation of hadronic current, 


^l(^) = Jl(j) = r+(gv.g^o + gASi,ia^j)E{q‘^) , (12) 

j 

where fj is the position of the j-th nucleon in the nucleus and E{q^) is the form factor 
dehned by 


nq") = (- 


(13) 
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with the value niA = O.SSGeV. By these integration, we obtain the energy conservation 
Ei = Ef + ei + 62 and the energy denominators. We hnd 
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where en = En — {Ei + Ef)/2 and 
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{Qw + ki + k2)fj,qwn' 

{Ew + en+ ^-^)Ew[{Ew + € 2 )^ - E^^Ew + ei + 62 )^ - E^] 

_l_ Qw^iiqw — ki — k2)fj.' 
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( 15 ) 

Hereafter, we consider the case where the chirality selection rule is satished, i.e., 
VlVr = 0 . So far we took (al) the S-wave function for electron wave function and 
thus the total angular momentum taken by electrons are 0 or 1 so that the 0 ^ —O’*" 
and 0^ —> transitions are allowed in general. Then, we used (a2) the non-relativistic 
approximation of the hadronic current because we are dealing with the allowed transition. 
Next, we make (a3) the closure approximation where is replaced by the average value 
< Sn >= jJ-oTne- Then, the sum of the intermediate states can be taken. Then, we hnd 

y < N,\jUr,)\n X n\Jl,(r,m > = Y.< Nl\JUn)\n >< n|4(o)|A'i > 

n n 

= < N,\Jl{rj)jl.{mNi > ( 16 ) 

SO that and enter in Sfi as their sum. 

Now, we concentrate on the 0^ —0^ transition. Since Jv{^j) is an parity even 
operator, only the mu = p' = 0 and fi = k, fj,' = k' {k, k' = 1, 2, 3) parts contribute to the 
O’*" —0+ transition because q is an odd parity operator. Now, we find that Cqo + T*oo 
and Ckk' + Dkk' are even functions with respect to the exchange of ei and 62 . Then, the 
hrst part and the (ei 62 ) part are combined as 

X y < N,\JUri)jUr,)\N. > (C„,, + £>„,) 

( 17 ) 

Now we use the identity Then, we hnd 
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where we used Ckk' + Dkk' = qkQk'iC + D). 

Now, we expand and with respect to the small quantities ei/E^, ei/Ew, 
~ ^w) ^he leading order terms. Here we assume > mw and used 

the fact that ~ a few MeV. Then, we obtain 
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Here we used the fact that the momentum g is effectively cut off by due to the form 
factor E{q^). 

Then, we use 
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After taking the spin snm and performing the phase-space integration, we hnd the 
half-life of the transition of the nentrinoless donble beta decay dne to the heavy composite 
nentrino for 0+ —> O’*" transition as 
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where Gqi is the phase space factor dehned in Eq. (3.5.17a) in the paper by Doi, Kotani 
and Takasngi[9]. The above expression is obtained by taking | i/l 1“^ + | |^= 1 which 

is valid becanse of the chirality conservation. 


3 Constraint from neutrinoless double beta decay 


In the following, we analyze the constraint on the conpling by nsing the experimental 
half-life limit of ^®Ge which is measnred by the Heidelberg-Moscow collaboration)^, 

T(0+ ^ 0+ : ^®Ge) > 5.6 ■ lO^V 90%c./.. (28) 


We derive the constraint on composite parameters from this data. By nsing the valnes of 
nnclear matrix elements obtained by Hirsch, Klapdor-Kleingrothans and Kovalenko[10], 

Mgt,n = 0.113, M'p = 3.06 ■ 10■^ M'gt = -7.70 ■ 10■^ M!p = -3.09 ■ 10■^ (29) 


the phase space factor Gqi = 6.4 ■ 10 /yr given in Ref.6 and Qa/Q v = 1-25, we hnd the 
half-life of the 0^ —O’*" transition for ^®Ge is 
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Then, we find the constraint 
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4 Discussions 


Firstly, we shall discuss about the excited neutrino mass. We consider the case that the 
composite scale is the same as the excited neutrino mass m*. In Fig.l, we show the 
limit of m* for Xw > 1- From this figure, we find 

m* > SAmw (Xw > !)• (32) 

Firstly, we compare our new result with the old result which is estimated by using 
the effective interaction in Eq.(3). By comparing Eq.(28) and Eq.(5), we find that the 
second term in the parenthesis is the same as the old one, while the first term is different 
by the cancellation among the poles due to u* and W propagator. 

Next, we compare our result with the one by Panella et al.[8]. One difference is 
the coefficients of mass factors of nuclear matrix elements. They used the assumption 
m* >> mw- Our coefficients agree with them in this limit. We evaluated without this 
assumption to see the behavior near mw- Our formula is valid for m* >> m^ = 0.85GeV. 
Another difference is the overall factor. Our decay rate formula is about sixteen times 
larger than their formula for m* >> mw case where we can compare with their formula. 
As a result, we find the factor 2 stringent limit on the composite scale m^^, or the relative 
coupling Xw- These can be seen from Fig.2 and Fig.3. In Fig.2, we showed the lower 
bound on Ac = mu^/^/2 for Xw > 1- The lower bound at m* = 6 mw is about 0.15 
while their bound is about 0.08. In Fig.3, we showed the upper bound on when 
Ac = m^^/\/2 =lTeV. Again our bound is about factor two severer than their bound. 

In summary, we computed the half-life formula of the neutrinoless double beta decay 
for the 0^ —> 0^ transition. We made the systematic analysis without assuming m* >> 
mw- Our result differs from the one by Panella et ah, some mass factors due to their 
assumption m* >> mw and the normalization difference about sixteen times. 
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Figure Captions: 


Fig.l; The lower bound on the excited neutrino mass m* for = m* and Xw > 1. The 
allowed mass range is the region m* > “iAmw- 

Fig.2; The lower bound on the composite scale Ac = my^/\/2 as a function of m^/mw- 

Fig.3: The upper bound on the relative coupling Aye for Ac = > ITeV as a 

function of m^/mw 
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